program testspeed

real(8)::a(1000,1000),I,TA(2)
integer::x,y,k
a=1D0
I = DTIME(TA)
write(*,*) TA

do k=1,1000
a=a/2D0+a/2D0
enddo
I = DTIME(TA)
write(*,*) TA

do k=1,1000
do x=1,1000
do y=1,1000
a(x,y)=a(x,y)/2D0+a(x,y)/2D0
enddo
enddo
enddo
I = DTIME(TA)
write(*,*) TA
 end

